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/iBSTRACT 

The critical discharge of two- component, two-phase mixtures 
through convergent nozzles is analytically investigated for in- 
itially dispersed or annular- dispersed flow patterns. The proc- 
esses of interphase heat and momentum transfer are considered, 
and based on the experimental results of previous investigators, 
credible assumptions are made for these processes. The result-’ 
ing model yields predictions for both the critical pressure 
ratio and flow rate as a function of the stagnation conditions. 
The analytical results exhibit good agreement with the published 
air-water nozzle data. 

The model is extended to compressible flows through orifices 
and nozzles and again exhibits good agreement with the available 
two- component data. 

INTRODUCTION 

The compressible flow of two-phase mixtures has received 
considerable attention because of its occurrence in such prob- 
lems as the blowdown of reactor pressure vessels [1], the vent- 
ing of rocket propellant tanks [2], the design of Jet pumps [3], 
and the metering and control of two-phase flow [4]. Generally, 
the practical problem involves the flow of a one -component mlx.- 
ture, however, much information about such. flows can be obtained 
from observing the behavior of two-component mixtures. For ex- 
ample, several experimental studies have employed air-water mix- 
tures to observe the interphase processes of heat and momentum 
transfer. This information is extremely difficult to extract 
from one- component flows because of the simultaneously occurring 
mass transfer process. 

The objective of this investigation is to closely examine 
the available air-water nozzle data and employ this information 
to construct a compressible flow model which incorporates rea- 
sonable approximations for the rates and amounts of interphase 
heat and momentum transfer. The result will be a useful guide 
in structuring a mathematical model for one- component, two- 
phase compressible flows. 

ANALYSIS 

The steady-state, one-dlmehslonal, continuity and momentum 
equations for two-component (x » constant) flows can be written 
as 

Liquid Continuity 

W^Vj« A2 _u^ (l) 
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Gas Continuity 


Momentum 


W^v_ = 
g g g g 


- A 


dZ 


= W 



+ Wi 


dui 

dZ 


+ 4 . d£ 

dZ dZ 


( 2 ) 


(3) 


Por the high velocities characteristic of critical flow, the 
frictional and hydrostatic head loss quantities can be neglected 
with respect to the inertial and pressure gradient terms. For 
simplicity, only isolated systems are considered. 


Equation 3 can be approximated by 


dZ 



+ (1 - X) 


duj“ 
dZ _ 


(4) 


It is assumed that, for fixed stagnation conditions, Ug and u^ 
are of the form f[P(Z)]. Hence, Eq. 4 shows 



^ + (1 - X.) 



(5) 


(The subscript t indicates that all the enclosed quantities are 
evaluated at the throat.) The critical flow criterion states 
that for fixed stagnation conditions, the flow rate exhibits a 
maximum with respect to the throat pressure. 


If k is defined by 


k - s= X 1 - cx g 
u^ 1 - X a vi 



the critical criterion can be applied to Eq.. 5 to give: 



rflr 




dv. 





1 )] 

( 8 ) 


For all but the extremely small qualities (x of order 10"^), the 
role of the liquid compressibility can be neglected, hence, the 
liquid phase is considered to be Incompressible. 


The terra (dVg/dP)^ represents the cbmpresslbility of the 
gaseous phase and may be expressed as (dVg/dZ)^/(dP/dZ)^. If It 
is . assumed that the gas is perfect^* (PVg = RTg). then. 


\ 


(dVg/dZ)^. / Vg 

“(dP/dZ)^ ~ ' 1 ^ 

Therefore, to describe the compressibility at the thi r>at one 
must describe the axial temperature gradient which is indicative 
of the Interphase heat transfer rate . Smith, et al. [5] meas- 
ured axial temperature profiles of a liquid film in the 
dlspersed-annular critical flow of air-water mixtures through an 
annular venturi. These profiles are shown in Flgv.l. . 



Smith purposely generated a liquid film on the straight 
wall (Fig. l) upstream of the convergent region. Photographs 
show that this wall film persists through the convergent and 
throat regions and into the divergent region. The schlleren 
photograph of compressible gas-particle flows presented by 
Gilbert, et al. [6] also suggest that the more dense phase is 
concentrated near the centerline in the throat vicinity. (For 
Smith’s geometry, the centerline would be equivalent to the 
straight wall.) Such preferential distributions are in agree- 
ment with what one would expect from the two-dimensional flow 
field in the throat vicinity. Based on these observations, a 
flow pattern development like that shown in Fig. 2 is assumed to 
be typical of convergent nozzle flows, and, thus. Smith’s results 
can be considered to be representative of all such flows. This 
flow pattern structure is only proposed for flows in which the 
inlet pattern is annular, annular-dispersed, or mlrt and not for 
bubbly mixtures like those observed by Muir and Elchhorn [7]. 


Table I compares the measured stagnation to throat temper- 
ature drops shown in Fig. l to the predictions of the thermal 
equilibrium model which assumes the gas and liquid temperatures 
are equal at all times. If the gas expansion is represented by 

a'poiytropic process (Pv = constant), the thermal equilibrium 
assumption gives * 


(l - x)cij, + xc 
■" (l - x)c^ + xc 


£ 

V 


( 10 ) 


This is the same formulation derived by Tangren, et al.[8]. 

Table I also includes measurements reported by Vogrln [9] for 
temperature differences between the two stations shown in Fig. 3. 
The accompanying radial void fraction profiles reported in ref- 
erence 9 indicate that the downstream measurement is representa- 
tive of a liquid wall film. Therefore, the measurements are in- 
terpreted in the same manner as those of reference 5. The data 
given in Table I and Fig. 1 show the liquid temperature Is es- 
sentially constant throughout these rapid expansions. Thus, 
there is no significant amount of heat exchanged with the gas- 
eous phase which, therefore, expands in an adiabatic manner. 
However, the temperature profiles of Fig. l show a large temper- 
ature gradient at the throat, hence, the rate of change is quite 
different from an adiabatic process. If it is assumed that, for 
a given phase distribution, the rate of interphase heat transfer 
is proportional to ,the temperature difference between the phases: 




N 


q - hAg(Ti “ Tg) (11) 

The amount of heat transferred in the expansion is the rate in- 
tegrated over the time necessary for the mixture to travel from 
the stagnation region to the throat. 



( 12 ) 


Figure 2 shows the one-dlmenslonal pressure profile for a typical 
nozzle configuration. In the convergent portion, the accelera- 
tion and the accompanying steep pressure gradients occur between 
the upstream location which has a diameter twice that of the 
throat and the throat Itself. The temperature profiles, which 
are based on the results of reference 5, show that where times 
are long (low velocity), the temperature difference is small 
and that where the temperature difference is appreciable, the 
residence time is short (high velocity). Therefore, it is as- 
sumed that the expansion is so rapid that the amount of heat 
transferred is negligible (Q « O). 

(13) 

Figure 2 also shows that the rate of heat transfer at the throat 
can be quite large. If It is assumed that the gaseous behavior 
in the throat region can be represented by emother polytropic 

process (PVg^^ = constant), the exponent n^ must reflect a 
sizable heat transfer rate and, thus, cannot be equal to the 
Isentropic exponent y. As a compromise between simplicity and 
the real process, it is assumed that n^ is equal to the equi- 
librium exponent (ng) given in Eq. 10, which is indicative of a 
large rate of Interphase heat transfer. 




(14) 


The above arguments have essentially separated the expansion 
into two parts such that one approximation can be made for the 
overall expane" and another for the rate process at the throat. 


To eval ae critical flow rate, it is necessary to de- 

scribe the derxvieive dk/dP)^; which is representative of the 
Interphase momentum transfer rate. Using the axial void fraction 
measurements of Vbgrin and the approximation given in Eq. 13, 
one can calculate the velocity ratio profiles shown in Fig. 3. 
Since the velocity ratio is assumed to be of the form k[P(z)], 


_ dk/dz)t 
dP/t ■" dfP/dzjt 


(15) 


Figure 3 shows dk/dz)^ is either zero (a minimum) or small 
(k nearly constant) while dP/dz)t is quite large. As a first 
approximation it is assumed that 



= 0 


(IS) 


dp/t 

for nozzle flows operating in the annular, annular-dlspersed, 
or mist flow ranges. 

The approximations given in Eqs. 14 and 16 and the incom- 
pressible liquid assumption (dv^/dP = 0) enable one to simplify 
Eq. 8 to: 


kn„P 


G? = f. "“E 


c ) U + x(lc - l) ixVgJ 


(17) 


If the velocity ratio is assumed constant throughout the 
expansion, the two-phase frictionless momentum equation (Eq. 4) 
can be written as 


k[(l - x)v^ + xvg.] 
i + x(k - l) 


dP = d(^ 


(18) 


Using the assumption given in Eq. 13, along with Eq. 17, this 
expression can be Integrated between the stagnation and throat 
locations to give: 


(1 - x)kvj(Po - Pt) +7^ (PoVgo - Pt^gt) 


-{ 


[(1 - x)kvj + xVg] 


t 


(19) 


Equation 19 can be rearranged and expressed more compactly as: 


T] 


j- - °t 


T) 


,i/y 


'"E + (1 - 'H) 
_aof “t 


+ 


7-1 


( 20 ) 


where 


*0 


( 21 ) 


The constant velocity ratio assumption is obviously not in agree- 
ment with the axial, profiles shown in Pig. 3. However, the high 
upstream velocity ratios shown occur in the stagnation region 
where system velocities are negligible and the velocity ratio is 
essentially meaningless. The assumption Implies the velocity 
ratio is essentially constant in the region of high acceleration. 
A similar assumption has been widely used for particulate com- 
pressible flows [10]. 


X 


Since the constant velocity ratio assumption dictates that 


°° 1 - at(l - 


(22) 


the critical pressure ratio is only a function of either the 
stagnation or throat void fraction and the poly tropic exponents 
y and ng. Figure 4 compares the proposed model and the ther- 
mal equilibrium solution with the experimental results of ref. 9 
as a function of stagnation and throat void fractions. It is 
seen that there is little difference between the two solutions, 
however, of the two models, the one proposed herein appears to 
be more characteristic of the data. 


Unfortunately, the void fraction is usually unknown, and, 
thus, it is preferable to express the critical pressure ratio 
as a function of quality. As shown by Eq. 19, such an expres- 
sion requires a knowledge of the throat velocity ratio. These 
velocity ratios can be calculated for Vogrln’s experiment by 
using the measured throat void fractions and the adiabatic ex- 
pansion assumption discussed above. The values, as shown in 
Pig. 5, appear to exhibit some pressure dependency, however, 
the uncertainty is so large that any attempt to discern such a 
relationship would be highly speculative. Therefore, as a first 
approximation, it is assumed that the velocity ratio is dependent 
upon neither the pressure nrr the quality and equal to the aver- 
age of the data shown (k-^^ = 3.2). Equation 20 can be written as 


T1 - 


(1 - x)ktVl _y_ 1/y 

xv-4, 7 - 1 ' 



r(l - x)ktvi ■ 

2 

+ 

"(1 - x)k-{jvf 

, y 

2 

L J 

L J 

' 7 - 1 


where 


RT 


gt 


(n-V-V) 


and the critical flow relation is given b: 


(23) 


(24) 


C II + x{kt - 1 )]xRTq 


(25) 


For given conditions of x, Tq, and either the stagnation 
or throat pressure, the transcendental critical pressure ratio 
expression can be solved. Q?hls solution implicitly Involves 
the critical flow rate, thus, a solution of Eq. 23 yields values 
for the critical pressure ratio and mass flow rate-. 


In summary a model is proposed for the critical flow of a 
two- component mixture which is initially in an annular, annular- 
dlspersed, or mist flow configuration. The model varies from 


X 


those existing in the literature because it employs separate 
assumptions for the amounts of heat and momentum transferred 
in the expansion and the rates of these Interphase processes 
occurring at the throat. The proposed one-dimensional solu- 
tion is dependent on the throat velocity ratio and is somewhat 
limited by the lack of experimental measurements of this 
quantity. 

COMPARISON WITH EXPERIMENTAL RESULTS 

It is convenient to nondlmensionallze the critical flow 
rate with the prediction of the well documented Homogeneous 
•Thermal Equilibrium Model. This model is presented in ref. 8 
and is based on the following assumptions: 

1. The average velocities of the phases are equal (k = l) 
at all times. 

2. Thermal equilibrium exists between the phases (T^ = T, ) 

at all times. ' ^ 


3, The system entropy is constant. 


The critical pressure ratio and flow rate predictions for this 
model are given by: 


and 



iL.:,. 

XVgH "E - 1 

n^ pi - x)v; ^ r(l - ^)vi' 

2 L ^gH J L *'"gH . 


+ 




xR^ 

o 



(26) 


(27) 


where 



'28) 


Dividing Bq. 25 by Eq. 27 gives 


^c 


T 

^cH " 

|^[1 + x(k|. - l)] 

■^"E+U/nE- 

am «■ 

[f 


(29) 


This egression is only slightly dependent on pressure and, thus, 
is a convenient means with which to compare the proposed model 
and experimental data over a wide range of stagnation pressures. 


The critical pressure ratio predictions of Eqs. 23 and 26 
are compared to the experimental results of Vogrin and Smith in 


Pig. 6. (For all the data and models presented herein, the gas- 
eous phase is assumed to he at 100^ relative humidity in the 
stagnation region and the weight fraction of water vapor remains 
constant at this value throughout the expansion. ) Here again it 
is seen that there is little difference between the Homogeneous 
Thermal Equilibrium Model and the proposed solution throughout 
the quality range considered. However, the data in the low 
quality region of annular-dlspersed flows are generally better 
described by the proposed model. The discrepancies at the higher 
qualities are due to the oversimplification of the Interphase 
heat transfer rate at the throat. The prediction of Chisholm [11] 
is also included to Illustrate the sensitivity of these relation- 
ships to the values chosen for the polytropic exponents. The 
solution presented by Chisholm differs from the proposed model 
in two ways. First, the velocity ratio correlation is formed 
from the long constant area duct results of Pauske [12], and 
secondly, the total expansion between the stagnation and throat 
regions is assumed to be isothermal and the throat polytropic 
exponent (n^) is taken to be 1.2. These thermodynamic approx- 
imations are essentially the inverse of those made in this study. 
The critical pressure ratio results given in Fig. 6 are mainly a 
function of the values chosen for the polytropic exponents, and 
the available data support the assumptions proposed herein as 
opposed to those of Chisholm. 

Figure 7 compares the normalized flow rate predictions of 
Eq. 29 and those of Chisholm with the data of Smith and Vogrln. 

For a stagnation pressure of 30 psia, the predictions of the two 
models are quite close, however, the Chisholm model exhibits a 
strong pressure dependency which is not characteristic of the 
data. Tne slight pressure dependency given by Eq. 29 is more 
representative of the experimental results.. 

By considering Figs. 6 and 7 together, the reader can dis- 
cern the merits of the proposed model. The Homogeneous Thermal 
Equilibrium Model yields a good prediction for the critical 
pressure ratio but it is off by as much as 100^ In the flow rate 
prediction. On the other hand, Chisholm’s approach gives a good 
estimate of the flow rate but an erroneous value for the pressure 
ratio. Therefore, neither of these two solutions correctly de- 
scribes the physical phenomenon. The proposed model yields good 
predictions for both the critical pressure ratio and flow rate 
and, thus, is more characteristic of the actual behavior. 

The geometries of refs. 5 and 9 where similar in that they 
were both straight wall converging-diverging nozzles. Smith re- 
ported results for flows with and without a liquid wall film up- 
stre^ of the convergent region. The data without a wall film 
apecclosfirttso the upstream flow configuration of Vogrin’s data 
and generally exhibits better agreement with the results of 
ref. 9. It is expected that the throat velocity ratio is de- 
termined by the flow pattern development which is dependent on 
the nozzle geometry and the inlet flow configuration. The data 
shown in Pigs. 6- and 7 demonstrate that the inlet flow pattern 
exerts a very small Influence on the system compressibility. In 
contrast, the data of Graham [13 j, shown in Fig. 8, clearly il- 
lustrate the influence of nozzle geometry on the critical flow 



rates. (Since Graham exhausted directly to the atmosphere, only 
runs with stagnation pressures greater than 30 psia are con- 
sidered. This ensures the results are characteristic of choked 
flow,) Since each of the geometries reported shows ^q/^cEE 
be essentially Independent of pressure, it seems reasonable to 
conclude that the discrepancies between Eq. 29 and the data are 
a result of the crude approximation given in Pig, 5. When the 
drastic differences between the test nozzles are considered, it 
is not surprising that this approximation is not universally 
valid. The model generally overpredicts the experimental re- 
sults, however, it does correctly describe the trends with qual- 
ity and pressure. This gives support to the assumptions employed 
to represent the amount of heat transfer and the rates of inter- 
phase heat and mass transfer at the throat. These assumptions 
are further Justified by comparing the model wlth'.;the high pres- 
sure air-water data of Netzer [l4j. Obviously, at high pres- 
sures, the throat velocity ratio will be smaller because the 
density difference between the phases is decreased. Table II 
compares the data of ref . 14 with the proposed model for throat 
velocity ratios of 1.0, 1.3 (the approximate value generated by 
Netzer’s numerical analysis), and 3.2. The model shows better 
agreement for smaller velocity ratios, which is the expected 
trend, and yields a good prediction for the critical pressure 
ratio which is only slightly dependent on the throat velocity 
ratio. Therefore, it is seen that the proposed solution is 
properly sensitive to the throat velocity ratio, but, as dis- 
cussed above, it is limited by the lack of throat velocity ratio 
data as a function of stagnation pressure and nozzle geometry. 

For the pressure levels shown in Fig. 8, flows with qual- 
ities in the range 0.005 < x < 0.010 should begin to transition 
into bubbly flow patterns, and the model is not Intended for this 
type of flow configuration. 

EXTENSIONS OP THE ANALYTICAL MODEL 

In addition to the experimental results for convergent 
nozzles, attention has been given to the compressible flow of 
two-phase, two-component mixtures through orifices and cylin- 
drical nozzles. If the single phase compressible behavior in 
such geometries is known, the two-phase model developed above 
can be extended to these configurations. 

Perry [15] has experimentally demonstrated that the dis- 
charge coefficient for the compressible flow of air through a 
sharp edged orifice is approximately 0.84. Jobson [16] has 
analytically shown that this value is consistent with the in- 
compressible coefficient of 0.61. It is assumed here that the 
single phase compressible value is applicable to two-phase com- 
pressible flows. The discharge coefficient for flow from a 
stagnation condition is defined by 



This definition alters the critical pressure ratio expression 
(Eq. 23) to 


(1 - x:)ktv i T, l/y 

7 - 1 ' 

i tr -* “^2 r ( 1 . v^lr. 


ng 

-(1 - x)ktvi , - 

2 

+ 

"(l - x)k^vf 

+ Y . 

2 q 2 

XV. 

L 6^ J 

L ^gt J 

■ Y - 1 


( 31 ) 


This relationship was solved for given conditions of x, Tq, 
and Pq with kt; = 3.2 as shown in Fig. 5. The resulting 
critical flow rate prediction is compared to the sharp ec ^ed 
orifice data of Graham [13] in Fig. 9. The good agreeme be- 
tween the m,odel and the data is readily apparent. Here ain, 
as the flow regime becomes a bubbly mixture, the propose model 
is not expected to be representative of the data. 


Guzhov and Medvedev [17] have reported data for the com- 
pressible flow of air-water mixtures through cylindrical nozzles 
(L/D ~ 4 and 5). The experimental discharge coefficients for 
single phase compressible flow through such configurations ap- 
pears to vary -from 0.80 to 0.93 [18:]. Therefore, as a first 
approximation, it is assumed that the orifice coefficient, 0.84, 
is also applicable here. The experimental data presented in 
[17] are for throat pressures ranging from 17.4 to 23.2 psia. 

For a basis of comparison, an average throat pressure of 20 psia 
was chosen for computation. Figures 10 and 11 compare the re- 
sults of Guzhov and Medvedev with the model proposed in this 
study. The model describes the data trends very well and ex- 
hibits good agreement for x > 0.004. This quality is in or 
near the region of a transition to a bubbly configuration, thus, 
it is not ■Burprl-3-ing-4^hat-t4ie-.jnadel_l^^ not representative of the 
data for qualities less than 0.004. ^ 


SUMMARY AND- CONCLUSIONS 


A model is developed for the critical discharge of two- 
component, two-phase dispersed, annular-dlspersed, or mist flow 
mixtures through convergent nozzles. The basic approach em- 
phasizes the need for separate approximations to represent the 
amount of interphase heat and momentum transfer throughout the 
expansion and the rates of transfer at the throat. 

The proposed solution is superior to those presently in the 
lltei“ature because it yields accurate predictions for both the 
critical pressure ratio and the critical flow rate. However, 
the general application of the model is limited by the lack of 
throat velocity ratio data as a function of the pressure level 
and no'zzle geometry. 

Based on the results of single phase flow, the model can be 
extended to the compressible flow of two- component mixtures 
through orifices and cylindrical nozzles. A comparison with the 
available data shows good agreement with the experimental results 
for such configurations. 


NOMENCLATURE 


n 



A area 

C discharge coefficient 

c specific heat 

Fw wall Shear 

f function 

G flow rate/unit area 

h heat transfer coefficient 

k velocity ratio Ug/u^-/ 

n polytropic exponent 

P pressure 

Q heat 

q heat per unit time 

R universal gas constant 

T temperature 

t time 

u velocity 

V specific volume 

W flow rate 

X quality Wg/W^ 

Z axial length 

a void fraction Ag/(Ag + Ai) 

P volume fraction 

y Isentropic exponent 

t] critical pressure ratio Pt/^o 

0 hydrostatic head force 

p density 

T transit time between stagnation and throat locations 

Subscripts 

c critical flow 

E thermal equilibrium (Tg = T^) 

g gaseous phase 

H homogeneous 

1 liquid phase 

m mixture' 

o stagnation 

p constant pressure 

B surface 

t throat 

V constant volume 
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TABLE I. - MEASURED LIQUID TEMPERATURE DROPS IN 

RAPID EXPANSIONS 


Investigator 

Quality, 

X 

AT (°P) 
Measured 

AT (Op) 
Thermal Equl. 

Vogrln 

0.0587 

0.0 

-3.42 


0.0684 

+0.5 

-4.45 


0.0817 

-1.5 

-5.76 


0.0891 

-0.5 

-6.76. , 

Smith, et al. 

0.945 

-1.5 

-55.5 


0.827 

- 2.0 

-45.5 


TABLE II. - COMPARISON OF ANALYTICAL MODELS WITH THE 
EXPERIMJENTAL DATA OP NETZER 


Experimental Data 
Horn. Ther. Equl. Model 

Proposed Model 

Chisholm Model 


Experimental Data 
Horn. Ther. Equi. Model 

Proposed Model 

Chisholm Model 


(ps?a) 

515 
515 
r 515 
< 515 
1515 
515 


515 
515 
f515 
< 515 
L515 
515 


0.091 

0.091 

0.091 

0.091 


0.130 

0.130 

0.130 

0.130 

0.130 

0.130 


1.0 

1.0 

1.3 

3.2 


1.3 

3.2 


T] 



(ibm/sec-ft^) 

0.55 

5730 

0.56 

4770 

0.57 

5248 

0.55 

5770 

0 48 

7482 

C.49 

4885 

0.57 

5050 

0.57 

4220 

0.58 

4504 

0.57 

4957 

0.52 

6410 

0.50 

4190 





















Figure 1. - Axia! liquid fitffl temperature profiles 
reported in Ret. 5. 
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Figure I - Proposed flow pattern deveiopment and accompanying 
lieat transfer beliavior. 
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Figure 5. - ThroU vehKity ratios interpreted from the ex- 
perimental results of Ref. P. 
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Figure 6. - Comparison of the critical pressure ratio 
predictions and the experimental results of Refs. 5 and 9. 


NORMAUZEO CRITICAL FLOW RATE.. NORMAUZED CRITICAL FLOW RATE, 

GrIGrUC OrlGrUC 


DATA 

A SMITH (UPSTREAM WALL FILM) 

26PSIA<P-s50PSIA 

* SMITH (NO UPSTREAM WALL FILM) 

26PSIAlPo<38PSIA 
o VOGRIN3OPSIA<P{,<80PSIA 

• VOGRIN50?SIA^P„<80PSIA 

® VO(,RIN80PSIAiPo<105PSIA 

Mooas 

PROPOSED MOOa 


CHISHOLM 



Figure 7. - A comparison of thecriticai iiow predictions 
and the experimental results at Refs. S and 9. 
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figure 8. - A comparison of the critical How rate prediction 
and the experimental results of Ref. 13. 




Figure 9. ' A comparison Mween the orifice critical (low 
prediction and the experimental results of Ref. 13. 



Figure 10. ■ A comparison at the critkal pressure ratio 
prediction and the cylindrical nozzle dMa of Ref. 17. 



Figure U. ~ A comparison between the critical (low pre- 
diction and the cylindricat nozzle data of Ref. 17. 


